CONNECTION BETWEEN THE SLAVE-PARTICLES AND 
X-OPERATORS PATH-INTEGRAL REPRESENTATIONS. A 
NEW PERTURBATIVE APPROACH. 

A.Foussats*, A. Greco*, C.Repetto*, O.P.Zandron* and O.S.Zandron* 

* Miembros del Consejo Nacional de Investigaciones Cientfficas y Tecnicas - Argentina. 

Facultad de Ciencias Exactas Ingenieria y Agrimensura de la UNR, 
Av. Pellegrini 250 - 2000 Rosario - Argentina. 

Abstract 



Abstract 

In the present work it is shown that the family of first-order Lagrangians 
for the t-J model and the corresponding correlation generating functional 
previously found can be exactly mapped into the slave-fermion decoupled 
representation. Next, by means of the Faddeev-Jackiw symplectic method, a 
different family of Lagrangians is constructed and it is shown how the cor- 
responding correlation generating functional can be mapped into the slave- 
boson representation. Finally, in order to define the propagation of fermion 
modes we discuss two alternative ways to treat the fermionic sector in the 

path-integral formalism for the t-J model. 
PACS: 75.10.Hk and 75.10.Jm 
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I. INTRODUCTION 



Accounting the relevance in describing the behavior of strongly correlated electron sys- 
tems, a renovate interest in the study of the supersymmetric generalisations of the Hubbard 
model appears in the last years. A complete review on strongly correlated electron systems 
together with its connection to high-T c is given in Ref. [1]. 

The Hubbard models based on the superalgebras spl(2,l), osp(2,2) or su(2,2) have been 
formulated by using several approaches^i. For instance, as it was suggestedEHl the super- 
algebra spl(2,l) could be useful to study the model in the limit of infinite on-site repulsion 
and with infinite-range hopping between all sites'. 

Many problems about correlated electron systems were treated in the framework of the 
decoupled slave-particle representations. Two of them are the most important: the slave- 
boson and the slave-fermion representations. The first one privileges the fermion dynamics, 
and therefore the slave-boson representation seems to be more adequated to describe a Fermi 
liquid stateffl. Instead, the slave-fermion representation seems to give a good response 
when the system is closed to the antiferromagnetic order^S 

An important question in order to understand the physics of the high-T c superconductors, 
is to solve how to go from one representation to the other. This is because in high-T c 
superconductors, both the Fermi liquid and the magnetic order states seem to be present. 

On the other hand, one of the main problems appearing in these kind of models is to 
define the dynamics of fermions in the constrained Hilbert space, when the double occupancy 
of lattice sites is excluded. In this case also a convenient representation is given in terms of 
slave-particlesEl. 

As it is known the slave-particle models exhibit a local gauge invariance which is de- 
stroyed in the mean field approximation. This local gauge invariance has associated a 
first-class constraint which is difficult to handle in the path-integral formalism. 

Another possibility to attack the problem was given in Refs. [16,17] by using generalised 
coherent states in the framework of the functional integral formalism. 

Recently, the t-J model was analysed in the context of the path-integral formalism'^. 
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Our starting point was the construction of a particular family of first-order constrained 
Lagrangians by using the Faddeev-Jackiw (FJ) symplectic method!!, in the supersymmetric 
version0'@. In this approach any decoupling is used, but the field variables are directly the 
Hubbard X-operators which verify the superalgebra spl(2,l). In this way, we always work 
with the real physical excitations. 

Next, by using path-integral techniques, the correlation generating functional and the 
effective Lagrangian were constructed. Moreover, we have proved that our path-integral 
representation can be directly related with that found in Ref.[23]. 

As mentioned above, one of the interesting and non completly solved problem present in 
this constrained system, is to study the fermionic sector when the double occupancy of lattice 
sites is excluded. In particular the role of the fermionic constraints, and so the fermionic 
dynamics in the constrained Hilbert space is a crucial problem that must be investigated. 

Therefore, one of the main purposes of the present paper is to go deep into the discussion 
of the different alternatives which allows us to define the fermionic propagator in the t-J 
model. 

In Ref.[19], in the framework of the perturbative formalism the Feynman rules with 
appropriate propagators and vertices were found. In particular, a discussion on the fermionic 
propagator was also given. 

In order to continue with the study of the fermionic propagation, different alternatives 
are exposed in the present paper. We show that there is another way that the one given in 
Ref. [19] to obtain fermionic propagation. This is done by working inside the path-integral 
and integrating out the two delta functions on the fermionic constraints. 

Moreover, other interesting point is to check our formalism with those obtained by means 
of the slave-particle representations. More precisely, we will show how our path-integral 
expression for the partition function (see equation (4.1) of Ref. [18]), written in terms of the 
Hubbard operators, can be mapped in the partition function coming from the slave-fermion 
representation. 

On the other hand, by following the FJ symplectic method, it is possible to show that 
a new family of first-order constrained Lagrangians written in terms of the Hubbard X- 
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operators exists. This family of classical Lagrangians is able to reproduce the Hubbard 
X-operators commutation rules verifying the graded algebra spl(2,l). As it can be shown 
this family of Lagrangians, totally constrained in the boson-like Hubbard V-operators, can 
be mapped into the slave-boson representation. 

The paper is organized as follow. In section II, the main results of Refs. [18,19] are col- 
lected. In section III, by analysing the change in the constraints structure of the t-J model, 
it is showed how starting from our path-integral expression for the partition function previ- 
uosly found, the partition function coming from the decoupled slave-fermion representation 
can be recovered. In section IV, by using the FJ symplectic method, a different family of 
first-order constrained Lagrangians is found. This family of Lagrangians corresponds to the 
situation in which the bosons are totally constrained. In such conditions it is possible to 
show how the correspondent partition function can be mapped in the partition function 
coming from the slave-boson representation. In section V, two alternative ways to define 
the fermion propagation are studied. In section VI, conclusions are given. 



II. PRELIMINARY AND DEFINITIONS 

In the t-J model the three possible states on a lattice site are | a > — | >, | + >, | — >. 
These states correspond respectively to an empty site, an occupied site with an electron of 
spin-up, or an occupied site with an electron of spin-down. Double occupancy is forbidden 
in the t-J model. In terms of these states the Hubbard V-operators are defined as 

Xf =| ia ><ip\ . (2.1) 

In Eq.(2.1), when one of the index is zero and the other different from zero, the corre- 
sponding V-operator is fermion-like, otherwise boson-like. 

The Hubbard V-operators verify the following graded commutation relations 

[Xf , X] 5 \± = Sij{S^Xf ± 5 a5 Xf) , (2.2) 

where the + sign must be used when both operators are fermion-like, otherwise it corresponds 
the — sign, and i, j denotes the site indices. 
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We assume that the family of classical constrained first-order Lagrangians in terms of 
the Hubbard X-operators can be written as follows 

L = a aP {X)X aP - VW . (2.3) 

In the FJ language^ the symplectic potential V^ ) is defined by 

V {0) = H(X) + \ a Vt a , (2.4) 

where A a are appropriate Lagrange multipliers, and so the constraints Q a are defined by 

#V(o) 

Therefore, the symplectic supermatrix associated to the Lagrangian (2.3) can be formally 
writtenil 

1 j , (2.6) 

where the compound indices A = {(«/?), a} and S = {(7<5), 6} run in the different ranges of 
the complete set of variables defining the extended configuration space, and | A | indicates 
the Fermi grading. 

Following Ref. [18], our starting point is to consider the following partition function for 
the t-J model written in terms of the four boson-like operators {X + ~ , X~ + . X ++ , X ) 
and the four fermion-like operators (X 0+ , X°~ , X +0 , X~°) 

Z = J VXf 5(n a ) 6(Q i2 ) 5(E i3 ) S(E tA ) (sdetM AB )f exp % J dt L(X, X) , (2.7) 

where L(X,X) is given by 

(1 + Pi)Ui - 1 



L(X, (2 _ Vi)2 _ Z U 2 ( X i +X t - X t~X- 

+ ~ E { X T X ? + x f x " a ) - » E *W - H t-A X ) • (2-8) 

where Ui = X^ + — X l and V{ = X^ + + X t . 

The Greek indices a,/3 takes the values {+, — ,0}, the index a takes the values {+, — }, 
and H t _j(X) is the usual t-J Hamiltonian 



H t _j = £ t {j XfXf + ~ £ JiJ XfX~T - \ £ ^ X r^f , (2-9) 

besides, in equation (2.8) a term depending on the chemical potential P was added. 

In Eq.(2.7) sdetM^B is the superdeterminant of the symplectic supermatrix M^b defined 
in (2.6), and the bosonic and fermionic constraints at each site i are given respectively by 

n a = X+ + + Xr- + Pi -1 = , (2.10a) 



Q l2 = X+-X-+ + ~ (X++ - Xr-) 2 - [1 - ~(X++ + Xr)] 2 + P* = , (2.10b) 



— i3 



X°+X+- - X»-X++ = , (2.10c) 



~ i4 = X+°Xr+ - Xf°X++ = . (2.10d) 

where was defined Pi = X° + Xf° + X^Xr . 

In equation (2.8), the Lagrangian coefficients as well as the constraints (2.10) were deter- 
mined by using the FJ symplectic method with the condition to reproduce at the classical 
level the generalised FJ brackets or graded Dirac brakets of the t-J model (see Ref.[18]). 

In particular the constraint (2.10a) deduced by consistency, is the completeness condition 
which must be verified by the Hubbard X-operators and plays an important physical role 
in this constrained model as it will be seen later on. At this stage it is important to remark 
that in the Dirac language^ the two bosonic constraints (2.10a,b) are second-class one. 

Now, it is useful to write the boson-like Hubbard X-operators in terms of the real com- 
ponents S a (a = 1,2,3) of a vector field S and the fermion-like Hubbard X-operators in 
terms of suitable component spinors (Grassmann variables)!^ 

X+ + = ±-(l- Pt )(s + S l3 ), (2.11a) 
*r = ^(1-Pi){s-S a ), (2.11b) 
X+- = j- s (l - Pi ){S a + iS i2 ) , (2.11c) 



Xr+ = ±- s (l- Pt )(S tl -iS l2 ) , (2.11d) 
XT = , X°~ = ** + , (2.11e) 

X+° = , X°+ = , (2.11f) 

where s is a constant and the hole density in the new variables writes pi = \l/* + \l/ i+ + \l/*_\E'i_. 
Accounting the fermionic constraints (2.10) it results (1 — Pi)(l + p%) = 1. 

The real vector field S can be identified with the spin only when p — 0, i.e., in the pure 
bosonic case. 

By using the second-class constraints (2.10a) and after the change of variables is made, 
the partition function takes the form 

Z = J VS a VS i2 VS i3 V^ ic V% c VX t Z>£ (sdetM AB )j exp (i J dt L eff ) . 

(2.12) 

where the quantity . is the super Jacobian of the transformation (2.11). 

The effective Lagrangian L e ff defined in Eq. (2.12), in terms of the new variables reads 

= I E SaS "-s* Sa + i e K*» + m E **** - H t -j 

+ E + Si + Sf 3 - s 2 ) + $ - iS a ) - *+(* + 

+ + - + Ss)) 6] , (2.13) 
where the Hamiltonian i/t_j of the t-J model writes 

Ht-j = E Uj*i*** ja + A E M 1 - - P.) [Si^i + ^2 + - s 2 ] . (2.14) 

*J> i,j 

In equation (2.13) the parameters Aj, 4* and ^ are respectively suitable bosonic and 
fermionic Lagrange multipliers. 

At this stage it is important to remark that our Lagrangian formalism is independent of 
the underlying lattice dimension. 

In the next section we are going to analyse the equation (2.7) in the framework of the 
decoupled slave-particle representations. 



III. SLAVE-PARTICLE REPRESENTATIONS. 



The standard way to construct the classical Hamiltonian formulation for slave-particle 
models and subsequently give the canonical quantisation is developed in Ref.[15] . The 
starting point is to consider the general classical first-order Lagrangian for n bosonic fields 
b a and m fermionic fields ft, defined on a lattice 

L(PI b a , fl f h ) = ~ J2(blb ia ~ b\ a b ia ) + \ J2(flf ib ~ flfib) ~ H{bl b a , fl f b ) . (3.1) 

i, a i : b 

Both, bosons and fermions fields, are submitted to the slave-particle first-class constraint 
at each lattice site i 

n* = E&U»+£y&/»-i = °- ( 3 - 2 ) 

a b 

Looking at the equations (3.1) and (3.2) it can be seen that when the index a takes the 
values ± and the index b takes only one value, the six fields (four boson and two fermion 
fields) define the slave-fermion representation. On the contrary, when the index a takes only 
one value and the index b takes the values ±, the six fields (two boson and four fermion 
fields) define the slave-boson representation. 

With the aim to confront our results with others previously obtained, in this section we 
consider the slave-fermion representation. In particular it is possible to confront the corre- 
lation generating functional (2.7) with that obtained from the slave-fermion representation, 
and this relation is not trivial. 

In our approach all the constraints are in the Dirac picture^ second-class, while in the 
slave-particle representations the constraint (3.2) is first-class. Thus, when the Hubbard 
X-operators are decoupled a local gauge symmetry is made evident. 

The starting point is the correlation generating functional (2.7) with the Lagrangian 
(2.8). 

By computing the (sdetMAB)f appearing in equation (2.7) we find 

(sdetM AB )j = -i^^- . (3.3) 

x°+x +0 

where pi evaluated on the constraints is written pi = . 

i 
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Integrating out the fields components X i , Xf and X { by using the delta functions 
on the constraints written as follows 



5(n i2 ) = 5{Xt + Xr - X+-X-+) = -^5(X r - X \+; ) , (3.4a) 



V0+ vH 

5(~ l3 ) = 5(X»+X+- - Xtxr) = Xt + 6(Xr - j^+T ) > (3.4b) 



5(E l4 ) = 5(X+°X-+ - Xf°Xr) = X^X- - ^^-) , (3.4c) 



and taking into account the equality (l + p i )8(X i ' + j ) = 5[(l — p i )(X i j + j )] = 

i i 

_n x +0 x~ + 

5(X i i ~xT^ — ) (coming from the property of the Grassmann variables), the partition 

i 

function (2.7) takes the form 

r X + ~ X~ + X 0+ X +0 

Z = j VXt + VXt VXr+ VX?> VX»+ 5{X+ + + ^ + - 1) 

x exp (i f dtL*(X,X)) , (3.5) 



where L* (X, X) is given by 



2Y^ 

+ + - H{X) . (3.6) 

i,a i 

As it is known the change of variables that allows to write the remaining five Hubbard 
X variables in terms of the fields variables in the decoupled slave-fermion representation is 
defined by 

Xt + = b\ + b i+ , (3.7a) 
Xf = bl + h- , (3.7b) 
X~ + = b\_b i+ , (3.7c) 
X? + = h + f} , (3.7d) 
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X+° = bin . (3.7e) 

^From equations (3.7) it can be seen that the five Hubbard X-fields are given in terms of 
the six fields of the slave-fermion representation, so it is necessary to introduce an aditional 
condition among the six fields of the slave-fermion representation to make possible the 
transformation. 

We assume the following general linear form for the conditions in each lattice site 

0* = E( G - fo » + G l 4) +Hifi- Hi f} + K l = 0, (3.8) 

a 

where Gi a , G ia , Ki are bosonic parameters and Hi , H\ are fermionic (Grassmannian) pa- 
rameters. 

As it was commented above, when the Hubbard X-operators are written in a decoupled 
representation a local gauge symmetry is made evident. Thus, from a constrained system 
with a set of second-class constraints, it changes into a constrained system with a first-class 
constraint, and therefore a gauge fixing condition must be imposed. Therefore, the equation 
(3.8) is not other than the gauge fixing condition which corresponds to the local gauge 
symmetry appearing in the decoupled representational. 

A convenient choice is to take in equation (3.8): Gi + = i, G\ + = —i and the remaining 
coefficients all zero, i.e equation (3.8) reads 

<f>i = t(b l+ - b\ + ) = . (3.9) 

Later on, in equation (3.5) is introduced as usual the unity 

1 = J Vbl Vb la Vf} Vft 5{Xt + - b\ + b l+ ) 5(X+- - b\ + b^) 8(Xr+ - b\_b l+ ) 
x 5(X° + - fjb i+ ) 5(X?° - blmfa) Ji , (3.10) 

where Jj is the super Jacobian of the transformation (3.7) and (3.9), and its value is J, = 

(h + + bl + ). 

By integrating out the five variables Xi, the partition function can be written 

Z = J Vbl Vb ic Vfi Vft S(Q,) 5^) Ji exp (i J dt L(b a , 6+ , / , (3.11) 

10 



It is easy to see that the super Jacobian Jj is equal to minus the determinant of the Dirac 
bracket constructed from the first-class constraint f2j and the gauge fixing condition (3.9), 
i.e —det[Qi , <P%]dj where the first-class constraint fij in the slave-fermion representation is 
given by 

n i = £ 6 U- + /i t /i- 1 = - (3.12) 

jo- 
in the equation (3.11) the Lagrangian L(b a , b\ , / , f') reads 

L(K , 6+ , / , /t) = i E( fo t^ - + ^ E(^A - fifi) - H <An ti> /«) • ( 3 - 13 ) 

i,<r i 

Therefore, the above considerations show that our correlation generating functional (2.7) 
can be mapped into the correlation generating functional (3.11) coming from the slave- 
fermion representational. This mapping is a consequence of the fermionic constraints present 
in our expression for the correlation generating functional (2.7). We can also conclude that 
it is not possible to relate our correlation generating functional (2.7) with that corresponding 
to the slave-boson representation. 

The next question is how to construct from the symplectic FJ formalism a new family 
of first-order Lagrangian by using the Hubbard X-operators of the graded algebra spl(2,l) 
as fields variables, in such a way that the results can be mapped in the slave-boson repre- 
sentation. The problem is solved in the next section 

IV. CLASSICAL LAGRANGIAN AND CONSTRAINTS. SLAVE-BOSON 

REPRESENTATION. 

By following Refs. [18,25], we assume that the family of classical first-order Lagrangians 
in terms of the Hubbard X-operators can be written as follows 

L = a aP {X)X^ - V<°> . (4.1) 

where the five Hubbard X-operators X aa ' and X 00 are boson-like and the four Hubbard 
X-operators X a0 and X 0a are fermion-like. In the present case the symplectic potential is 
V<°> = H(X). 
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The Lagrangian functional coefficients a a p{X) that a priori are unknown must be de- 
termined by consistency in such a way that the graded algebra (2.2) for the Hubbard X- 
operators is verified. By following the steps of Ref . [18] it is straightforward to construct the 
symplectic supermatrix associated to the Lagrangian (4.1). Thus the symplectic supermatrix 
Mab is written in the form 

Abb B b f 

. C fb D ff , 

The Bose-Bose parts A bb is a (10 x 10)-dimensional matrix and it takes the form 





M AB = 



(4.2) 



dago 



dam 



A bb = 



dX a °' dX°' ! 



I da oo> 

dX™' ^ 'dX w 



\ 



dX<™ 
_d£2oo 



8X° 



8X00 
dX 00 





9^oo 


QX. a " (T "' 


Q X o" a'" 




dfloo 


dX m 


ax 00 















(4.3) 



The Bose-Fermi parts B b f (the Fermi-Bose parts Cfb = —Bjf) is a (4 x 10)-dimensional 
rectangular supermatrix given by 



B 



bf 



9a aa 9 a a"a'" 


9a^ da^u^n 


Q X o"o"> dXO" 


Q X °"°'" dX°° 


da aa daoo 


da a o dam 


8X00 dX a ° 


dX°° dX"0 


dn a „ a ,„ 




8X0" 


dX° a 


dQ 00 


dQ o 


dX°° 


dX"" 



(4.4) 



The Fermi- Fermi parts Dff is the (4 x 4)-dimensional matrix given by 
Dff = 



da Qa 
dX°°' 


T Q X 0o 


da a{) 
dX a °' 


-r 9X ao \ 


daoa 
dX"'° 


1 da cr'Q 

dX°° 


da a o 


_i_ 9a CT , 1 
-r Q X a0 / 



(4.5) 



where Q aa > and Q 00 are the appropriate bosonic second-class constraints defining the struc- 
ture of the constrained model. 

Once the symplectic algorithm is applied and the correspondent differential equations 
are solved the solution we found is 



QiOtr — 



2X?° 1 ' 



a ia0 — 



h V-Ocr 

2X?° 1 



(4.6) 



and the boson-like Lagrangian coefficients are all zero. 
The set of bosonic second class constraints is given by 
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V°° = X™ + X+ + + Xr-l = 0, (4.7a) 



nf = xf - -^f- = • ( 4 - 7b ) 



In particular the contraint (4.7a) is the completeness condition necessary to avoid the 
double occupancy at each site. 

In these conditions the symplectic supermatrix is invertible and the matrix elements of 
it inverse gives the correct Hubbard graded brackets (2.2), i.e 

(M^)- 1 = -i(-l)\ £A \ [A , B] , (4.8) 

where | e A | is the Fermi grading of the field variable A. 

Consequently the dynamics in this condition is given by the Lagrangian 

L(X, X) = -^o E(^° CT X ? + X °° X t) ~ H{X) . (4.9) 

The Lagrangian (4.9) together with the bosonic constraints (4.7) correspond to a situ- 
ation in which the bosons are totally constrained and the dynamics is carried out only by 
the fermions. 

The partition function corresponding to this solution reads 

Z = j VXf 5[X™ + Xt + + X r - 1] 5[Xf - -^o^} (sdetM AB )f 

x exp (i [ dt L(X,X)) . (4.10) 



By computing the superdeterminant of the symplectic matrix appearing in (4.10) it 
results 



{sdetM AB )f = ( detA \det(D - CA^B)] ) 2 = (X™) 2 . (4.11) 



Now, in order to confront the correlation generating functional (4.10) with those coming 
from the slave-boson representation some algebraic manipulations are needed. 
The first step is to make the following change of variables 

<p il = X°°-btb i = 0, (4.12a) 
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Vla o = Xf-flb l = 0, (4.12b) 

^ 0CT = X?° - f ia b\ = . (4.12c) 

Analogously to what happens in the slave-fermion representation, in the decoupled slave- 
boson one also an additional condition among the fields is needed. From the general linear 
equation (3.8) we choose by simplicity the following reality condition 



b\ -bi = 0, (4.12d) 



The super Jacobian Jj of the transformation (4.12) is given by 



= _ K Ui-ru tJ ^ 



Now, by introducing in the equation (4.10) the unity 

1 = Jvfl Vf la Vb\ Vb t J, 8{Xf - flh) 5(X°° - f ia bl) 5{b\h - X?) 5(b\ - b t ) , (4.14) 

and integrating out all the fields variables Xf^ 3 , after same algebraic manipulations it is 
possible to show that the partition function (4.10) takes the form 

Z = J Vb\ Vb t VflVf* 6^) 8(fa) (h + b\) exp (z J dt L(6+ , b , f a , f a )) , (4.15) 

where and fa are respectively the first-class constraint and the gauge fixing condition in 
the radial eaue appearing in the partition function of the slave-boson representationS. 
They respectively read 

n i =6tfe i +x;/i/«r-i = o. (4.i6) 

fa = i(bi-bt)=0. (4.17) 

Again, we note that the factor (b + b') in the equation (4.15) is precisely the value of the 
det[Qi , fa]o appearing in the gauge theories containing first-class constraints. 
The Lagrangian L(tf , b , fl , f a ) defined in the equation (4.15) is given by 
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L(b\ b, fl f a ) = - - bhk) + - J2(flU - flfia) ~ H . (4.18) 

i i,a 

In summary, from our approach and working with the Hubbard X-operators without 
using any decoupling representation, the new family of Lagrangians (4.9) is obtained. The 
respective correlation generating functional (4.10) is mapped into the solution provided by 
the slave-boson representation. It is important to note that the new path- integral (4.10) in 
terms of the Hubbard X-operators up to our knowledge was developed in the present paper 
for the first time. 

We can see once more how a second-class constrained model written in terms of the 
Hubbard X-operators, when written in terms of the decoupled slave-particle representations, 
is transformed in a constrained system where a local gauge symmetry is made evident. 

We think that the results we can obtain by using the partition function (4.10) with the 
Lagrangian (4.9) can be useful for regimes where the system is close to a Fermi liquid state. 

In a forthcoming paper the partition function (4.10) is studied in detail within the con- 
text of the perturbative formalism. Having in mind the difficulty to treat the first-class 
constraint inside the path-integral (4.15), our purpose is to construct the Feynman rules 
and the diagrammatics starting from the path-integral (4.10). Once an appropriate fermion 
propagator can be found, our first objective will be to analyse the properties of the fermion 
spectral function. 

V. TWO ALTERNATIVE WAYS TO DEFINE THE PROPAGATOR OF 

FERMION MODES. 

In this section we discuss two alternative ways to treat the fermionic sector in order to 
define the propagation of the fermion modes. As commented above, a crucial problem in 
the t-J model is to define the fermionic propagation in the constrained Hilbert space, when 
the double occupancy is forbidden. 

With the purpose to study this problem, in Ref . [19] the correlation generating functional 
(2.12) was considered at finite temperature by means of the " Euclideanization procedure". 
Moreover it was assumed that we are close to an undoped regime where the system is an 
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antiferromagnetic insulator. Under this condition there is a small number of holes and it 
can be assumed that the hole density pi =< pi >= constant. The constant value p of the 
hole density must be determined later on by consistency, for a given value of the chemical 
potential p. 

In these conditions, it is possible to treat the non-polynomic Lagrangian (2.13) in the 
framework of the perturbative formalism, and so it can be partitioned as follows 

L eff = L B (S, A) + L F (rj) + Z/(S, 77) , (5.1) 

where 

L B (S, A) = -f E S " S *- S * S " + 2s ' £ \S i3 
is i s + s i 

+ ^2 X J' [^i<%+-0i ~~ Si2S( i+I ) 2 — Si 3 S( i+I ) 3 + S 2 , (5.2a) 

and 



L F (r7) + z/(s, r,) = E K** + »H + E 1****1, 

i,a i,<r i,j,a 

+ J2 f l i M i r 1i . (5.2b) 



Considering the bilinear bosonic part of the equation (2.13), for a constant value of the 
hole density and taking J^- = constant, we arrive at the equation (5.2a). In the equation 
(5.2a) was defined J' = J(l — p) 2 . Moreover, the symbol J2i indicates sum over nearest- 
neighbor sites. 

Besides, in equation. (5.2a) it was assumed that the vector S is written 

S = (0,0,s') + (£, St, S3) (5.3) 

where S± , S2 , S3 are the fluctuations. So, the equation (5.2a) corresponds to the lowest 
order of the system fluctuating around on antiferromagnetic state. Moreover, we must 
consider s 7^ s because as known the local magnetisation in an antiferromagnetic state 
reduced from its classical value, even for the pure Heisenberg model. The value of s must 
be determined also by consistency. 
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In equation (5.2b) the four-component spinor rj 



, is constructed from the two 



spinors St and The physical two-component spinor 



M = 



(5.4) 



is restricted by the 

fermionic constraint equations (2.10c,d), and the two-component spinor £ = | is a 

Majorana spinor. 

In the same equation the 4x4 dimensional matrix M. is defined by 

Is + S.<x 
Is + S.o- 
where cr = ( y <J 1 ,a 2 ,<T 3 ) are the Pauli matrices. 

In this regime the diagrammatics and the Feynman rules can be found. In particular 
the bilinear part of the bosonic sector written in equation (5.2a) gives rise to the usual 
antiferromagnetic magnon propagator (see Ref.[19]). 

The bilinear fermionic part L F (r)) of the equation (5.2b) can be written in terms of the 
four-component spinor rj and it is given by 



(5.5) 



where in the Fourier space the symmetric non-singular 4x4 dimensional matrix G^j is 
defined by 

f-(iu n + fx) e k l(s + s') \ 

e k -(iu n + ii) l(s-s') 



o 







\(*-<f) 



f 



9 



9 



-f J 



KVn,V n ) 



(5.6) 



In Eq.(5.6) the quantities k and v n are respectively the momentum and the Matsubara 
frequency of the fermionic field, and was defined e k = —t J2i exp(— il.k). 

The functions / and g appearing in Eq. (5.6) are totally arbitrary. As easily can be 
seen, these functions do not appear in the Lagrangian due to the Majorana condition on the 
two-component spinor $ > . 
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In this scenary, the symmetric matrix defining the fermionic free propagator G(o)«/3 is 
given by the inverse of the matrix (5.6). The physical components of the free propagator are 
given by the matrix elements G(o)ii,(j(o)i2,(j(o)2i and G(o)22 and they were explicitly written 
and analysed in section IV. of Ref. [19]. So, the Feynman rules propagators and vertices are 
given straightforward and therefore the boson and fermion self-energy can be computed. 

However, some particular features of the fermionic free propagator must be emphasised. 

The trick of introducing an auxiliary two-component Majorana spinor is a way to obtain 
a free functional G(o) that really propagates physical fermionic modes. 

The electron spectral function is defined from the fermionic propagator G(o) Q/ 3 by con- 
sidering the components G(o)n and G(o)22- The matrix elements directly connected with the 
electronic properties, as for example the Fermi surfase (FS), are precisely G(o)n and G(o)22- 
The electronic spectral function measured in photoemission experimentsi! must be related 
with minus the imaginary part of these matrix elements. 

Contrary to the fermion propagator obtained by means of the standard Green function 
methodS, our fermion propagator contains two poles. It is important to say that for a 
given filled factor the chemical potential we obtain is exactly the same obtained by using 
the standard Green function method. 

By ploting the electron spectral function (see Fig. 1 Ref. [19]) we can see that the peak 
at negative energy must be interpreted as the extraction of an electron, while the peak at 
positive energy represents the addition of an electron to the system. Therefore, the two 
peaks account for the photoemission and the inverse of the photoemission respectively. The 
presence of these two peaks implies that for a given value of k, the state is not completely 
filled or empty. Photoemission experiments are only sensitive to the first peak. Then the 
first peak of our propagator must be related with the excitation measures in photoemission. 

An alternative way to treat the fermionic sector is to start from the equation (2.7) and to 
work inside the path-integral. So, by integrating out the two delta functions on the fermionic 
constraints S i3 and S^, the partition function can be written as follows 

Z = J VXf VX° + VX+° 6(na) 5(n i2 ) (sdetM AB )f exp i J dt L*(X,X) , (5.7) 
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where L* (X, X) is given by 

+ {^ + Xr + X?°Xp) - H(X) . (5.8) 

The total Hamiltonian H is defined by 

H = H t _ J + ^J2 X^Xf, (5.9) 



where the Hamiltonian H t -j defined in equation (2.9) must be evaluated on the fermionic 
constraints E i3 and H i4 . 

Due to the non-linearity of the constraints (2.10c,d), when the path-integration on the 
two fermionic fields Xf~ and X^° is carried out, the non-polynomial structure of the kinetic 
fermionic part of the Lagrangian is made evident, as it can be seen from the equation (5.8). 

After the four boson-like X-Hubbard operators are related to the real components S a 
(a — 1, 2, 3) of a vector field S and the remaining two fermion-like X-Hubbard operators 
written in terms of suitable component spinor fields (see equations (2.11)), the correlation 
generating functional (5.7) takes the form 

Z = J VS a VS i2 VS i3 VV*_ V\ (sdetM AB )f (^pj exp (i J dt L eff ) . (5.10) 

Now, the Lagrangian L e ff defined in equation (5.10) is given by 

^ -sB'-a) (^#) + *-*•-) 

- H + Y.MSl + S^ + Sl-S). (5.11) 

i 

where in equation (5.11) and hereafter the tilde over the fluctuations is omitted. 

The Hamiltonian H written in terms of the real vector variable S and the two spinor 
component fields and reads 

H — ^2 7 i c V 1 O N ^il^j'l + ^i 2 Sj 2 + S i3 Sj3 + S 2 + s(S i3 + Sjs) 

i,3 ^ S ^i3)y S + ^JZ) 

+ i(S n S j2 - SaSji)] V-iVtj + ^ E J ij(l - P0(1 " Pi) [SnSji + S i2 S j2 

i,3 

+ s a s j3 - s 2 } + 2,^ (~^:) • ( 5 - 12 ) 



(3 
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Again, the path-integral (5.8) is considered in the framework of the perturbative formal- 
ism at finite temperature, and we assume that we are close to an undoped regime (antifer- 
romagnetic insulator). 

After a rotation of spins on the second sublattice by 180° about the Si axis is performed 
the Euclidean and rotated Lagrangian Lfjj is obtained, and so the lowest order of the 
effective Lagrangian (5.11) can be partitioned as follows 



Lf f R f = L B (S, A) + L F (^_ U J + L 7 (S, vT J , 



(5.13) 



where 



and 



L B (S, A) = -1(1 - p) £ S " S * + S* S " + 2s' £ Vfc 

+ g^2 E J' ~ S i2 S( i+ i) 2 — S i3 S(i +I } 3 + S- 



s + s i s + s' i 



(5.14a) 



(5.14b) 



L J = £ _^L_ _ iS . 2 + 5jl + iSj2 ) 

1,3 



(5.14c) 



At this stage it is important to note that our effective theory does not contain fermion 
dispersion. This feature is also present in the slave-fermion theories when the fermion 
dynamics is generated via de interaction with virtual magnons. 

By making a Fourier transformation it is possible to see that the bilinear bosonic part of 
the Lagrangian (5.14a) allows to recover the structure of the bosonic propagator (antiferro- 
magnetic magnons) given by 



( i5%(l-7.) 

2u> n 



(s+s')d (0) 








ii%( 1 +7«) 








1 

2s' 








\ 



1 

2s' 

J'z(l--y q ) , 

16sV 2 / 



S(u n ,u' n ) 



(5.15) 



In Eq.(5.15) the quantity <i( ) is given by 
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s + s' 



q 1 n J ' 



(5.16) 



where the frequency u 2 is defined by 



0J„ 



zf ( s + sf 



:i-7). 



(5.17) 



As 2 \2{l-p) i 

Moreover, in Eq.(5.15) z is the number of nearest neighbors and was defined the quantity 
lq = jE/exp(«I.q). 

Analogously, the bilinear fermionic part (5.14b) reads 



where we have named 



L F = £*-(MtO Go 1 *-(fc,fn), 

k,u n 



2s 



(5.18) 



(5.19) 



The inverse of this scalar function given by 

s + s' 1 

Lrn 



(5.20) 



2s zV„ — /i 

is a (non-propagating) functional which only depends on the Matsubara frequency v n . 
Finally, in the approximation that we consider the unique three-leg vertex is defined by 
/ e(k')+e(k) \ 



U a 



s + s' 



i{e{k') - e{k)) 




5(q + k — k') 5(uj + u- v') 



(5.21) 



At this point, the problem is to analyse the bilinear fermionic sector, in order to give 
the prescriptions for the propagation of the fermionic modes. The usual way to solve the 
propagation of fermions is by means of the Dyson equation. As known the Dyson theorem 
allows to compute the inverse of the corrected fermion propagator in terms of the free fermion 
propagator and the self-energy. Therefore the propagator G(k, v n ) = [Gq (u n ) — T,(k, z/ n )] _1 
can be calculated in a straightforward way within the self-consistent Born approximationili§. 
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By using standard techniques the following expression for the self-energy at zero temperature 
is found 

y(lc . , (i + p), 2 2 v Ni - (i - il) l/2 \ 1/2 - 7 fc+g [i + (i - ilf>^] 2 

z Y 

X ^77— : r (5-22) 

iv n - uj q - n - h(k + q , %v n - Ug) 

The expression (5.22) is useful in the strong coupling case (t > J). On the other hand the 

self-consistent solution of this equation is necessary in order to obtain fermionic propagation, 

and it must be performed numerically. Once an appropriate self-energy function iv n ) is 

found the propagator G(k,u) remains well defined and it is possible to compute numerically 

the spectral function defined by A{k, v) = —^lim £ ^ G(k,u + ie). Finally, as well known, 

the correction to the bosonic propagator is given by 



ub 



-1 



^o] ab -^a b ' , (5.23) 



where D^ ab is the inverse of the free bosonic propagator (5.15) and the bosonic self-energy 
U a b reads 



(5.24) 

As it can be seen from the above equation the dimension of the underlying lattice and the 
physics depend on the parameters z, ^ q and though our general formalism is dimensional 
independent. 

It is important to confront our results with other previously given in the literature related 
with the spin-polaron theories!!. Like in these theories our starting point was to assume an 
antiferromagnetic order state. This physical assumption is directly connected with the fact 
that at lowest order the fermion is not propagating (see equation (5.20)). Then, in order to 
describe a metallic phase where the holes move coherently on the lattice, it is necessary to 
solve the self-consistent equation (5.22). 

The solution of our equations (5.22) and (5.24) together with a quantitative comparison 
with the spin-polaron theories is an important matter that deserves further study. 
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Another point to take into account in a future work is to study the relationship between 
our matricial propagator G(o) a /3 ( see equation (4.13) of Ref[19]) and those obtained by solving 
self-consistently the equation (5.22). 

As it is known, any model or approach will be considered as a good candidate to describe 
high T c superconductors when it is able to answer the question of why the antiferromagnetic 
long range order disappears for small values of doping (for intance p = 0.04 — 0.05). In the 
last years this problem was attacked from differents approachesli!00. In a future work and 
from our formalism, we will also have a response to give about this important point related 
with the disappearance of the antiferromagnetism. 

In the present section, the magnetic excitations that we have considered are antiferro- 
magnetic magnons and in addition we have assumed a strong long range antiferromagnetic 
order. Therefore, our next step must be to study the instability of the antiferromagnetic 
order and to analyse against which phase this is unstable. In order to have some idea about 
this fact, in Ref.[19] we have studied the magnon self energy effects on the magnetic spec- 
tral function. Besides the sofltening of the antiferromagnetic magnon we have also found 
a reduction of the magnetic spectral signal. These results were got by using our two pole 
bare fermionic propagator. In the near future and in order to improve our calculation we 
will solve a self energy coupled problem for both magnetic and electronic dynamics. 

VI. CONCLUSIONS 

As showed first in Ref. [25] for the pure bosonic case (su(2) algebra), in a classical La- 
grangian formalism it is not possible to introduce the full Hubbard algebra by means of 
constraints. Consequently, in a path-integral formulation it cannot be introduced the com- 
plete information of the Hubbard algebra namely, the commutation rules, the complete- 
ness condition and the multplication rules for the Hubbard X-operators. So, the Heisen- 
berg model treated in the Lagrangian picture only admits two second-class constraints, and 
these are the completeness condition X ++ + X + X 00 = 1 and the non-linear constraint 
X+- X~+ + \{X ++ - X—f = s 2 . 
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The above last constraint in not really the group Casimir operator. It can be shown that 
the presence of such constraint is consistent with the quantisation of a spin system in the 
limit of large spin s, or equivalently for magnetic order state. 

A similar situation actually occurs in the case in which the Hubbard X-operators verify 
the graded algebra spl(2,l), but in this case at least two solutions are possible. When the 
Hubbard X-operators closes the graded algebra spl(2,l) the t-J model described in terms of 
a first-order Lagrangian has the following possible solutions: 

i) One is the family of first-order Lagrangians (2.8) together with the set of second-class 
constraints (2.10). Two of them are bosonics and the other two are fermionics. In particular 
the constraint (2.10a) is the completeness condition. As it was shown the corresponding 
path-integral formalism is mapped into the decoupled slave-fermion representation. So, in 
this case our correlation generating functional (2.7) privilege the magnon dynamics of the 
system with a strong feature of magnetic order state (consistent with the large s nature of 
the constraint (2.10b)). 

ii) The different family of first-order Lagrangians (4.9) together with the new set of 
second-class constraints (4.7) is also a possible solution. In this case all the constraints are 
bosonic, and (4.7a) is again the completeness condition. As it can be seen the remaining 
four constraints (4.7b) are related to the multiplication rules. In this situation the bosons 
are totally constrained and the dynamics is carried out only by the fermions. Moreover, we 
note that not any non-linear constraint of the type (2.10b) appears. As it was shown the 
path-integral formalism corresponding to this dynamical situation is mapped in the slave- 
boson representation. Therefore, in our correlation generating functional (4.10) the fermion 
dynamics with a strong feature of Fermi liquid is priveleged. 

It is possible to conclude that once the set of second class-constraints is chosen, different 
families of Lagrangians are obtained, and so we can insure that each family contains different 
physics. 

Worthwhile to remark that both Lagrangian formalism are independent of the dimension 
of the underlying lattice. 

Moreover, as it was seen in all the cases the completeness condition appears as necessary. 
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As it is well known the completeness condition involve an important physical meaning. Such 
condition avoids at the quantum level the double occupancy at each lattice site. 

Finally, in section V two alternative way to define the fermion propagator were developed. 
By means of the trick of introducing auxiliary Majorana spinors, a free fermion matricial 
propagator having two poles was found. Later on, by integrating out two of the fermions 
using the delta functions, it was possible to obtain the non-propagating scalar function 
(5.20) in the fermionic sector. ^From this free "propagator" and by means of the Dyson 
equation the fermion self-energy can be evaluated straightforward within the self-consistent 
Born approximation. 

Since our path-integral formalism is mapped in the slave-fermion formalism and taking 
into account that the equation (5.22) for the self-energy at zero temperature has a similar 
structure to that obtained from the spin-polaron theories^, both results allows us to insure 
that our approach is consistent. 
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